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f^ ', We have analyzed the X-ray spectra of the microquasar GRS 1915+105, as observed 

D I with the PCA on the Rossi XTE, during periods of stable weak emission, outbursts and 

rapid flaring. We find that the complicated X-ray intensity curve of this source can 
be described by the rapid removal and replenishment of matter forming the inner part 
of an optically thick accretion disk, probably caused by a thermal-viscous instability 
^ I analogous to that operating in dwarf novae, but here driven by the Lightman-Eardley 

instability. We find that the mass accretion rate in quiescence is ~ 10~^ Mq yr~^. Only 
O ' ^ small fraction of the energy liberated by accretion is emitted as radiation. We suggest 

^ ! that most of this energy is advected into the black hole in the high-viscosity state of the 

j^^ ' outburst cycle. 
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1. INTRODUCTION 

The X-ray source GRS 1915+105 was the first 
galactic object to show superluminal expansion in ra- 
dio observations ( Mirabel fc Rodriguez 1994 ). The 
standard interpret ation of thi s phenomenon in terms 
of relativistic jets ( Rees 1966 ) unambiguously places 
the source at a distance D — 12.5 kpc, wit h the jet 
axis at an angle i ^ 70° to the line of sight ( Mirabel 



fc Rodriguez 1994]) . GRS 1 915+105 was discovered as 
an X-ray tr ansient in 1992 ( [Castro-Tirado, Brandt, fc 
Lund 1992 ). It was found to be in a bright state with 
the Rossi X-ray Timing Explorer (RXTE) in April 
1996, and has been bright in X rays since. During 
this period it has shown a remarkable richness in its 
variability behavior. It showed quasi-periodic oscil- 
lations, with frequencies ranging from 10^'^ Hz to 67 
Hz ( Morgan, Remillard, fc Greiner 199^ ). On several 
occasions its light curve showed a complicated pat- 
tern of dips and rapid transitions between high and 
low intensity, which repeated on a time scale between 



30 and 400 minutes (Greiner, Morgan fc Remillard 
1996). The source is suspected to be a black hole 
binary because of its similarities to the other galac- 



tic superluminal source GRO J1655-40 (Zhang et al 



1994), which has a dynamical mass estimate imply- 
ing a black hole ( Bailyn et al. 1995 ), and because it 



has often been well above the Eddington luminosity 
for a neutron star. Here we show that the compli- 
cated light curve of GRS 1915+105 can be described 
by the rapid (~1 s) appearance and disappearance of 
emission from an optically thick inner accretion disk, 
probably caused by a thermal- viscous instability anal- 
ogous to that operating in dwarf novae. 

2. X-RAY OBSERVATIONS 

The observations discussed here were carried out 
with the PGA on the Rossi XTE, on 1996 June 19, 
September 27 and October 1 to 29. During these ob- 
servations the source was both very luminous (Lx — 
5 X 10^8 - 3 X 10^9 erg s'^ at D = 12.5 kpc in the 
2-50 keV band) and very variable (see Fig. 1), similar 
to the observations reported by Greiner, Morgan, fc 
Remillard (1996) and Morgan, Remillard, fc Greiner 
(1997). 

The variability is highly complex, but quite struc- 
tured. We identify three types of variability ( "states" ) 
from the fight curves (see Fig. 1): a quiescent state at 
a relatively low count rate and with little variability, 
an outburst state at high count rate with strong red- 



noise variability, and a flare state where the flux shows 
rapid alternations between two flux levels which can 
be extremely complex in shape. When more than 
one state is observed within the same observation (a 
typical observation spans a few hours), the sequence 
is always the same: quiescent-outburst-flare-quiescent 
etc., sometimes with a highly repetitive pattern. In 
some observations, the source is observed in only the 
quiescent or flaring state, but never only in the out- 
burst state. These states are similar to the lull, flare 
and sputter state described in Greiner, Morgan fc 
Remillard (1996) on the basis of previous observations 
of the source. 

In this Letter we take the observation of 1996 Octo- 
ber 7th (6:36 to 7:36 UT; Fig. 1) as a prototype con- 
taining all tfiree states. Tfie highly-organized com- 
plex variability also shows up in a plot of hardness 
ratio versus time (Fig. 1, bottom panel). The count 
rate and hardness ratio histories can be combined in 
a hardness-intensity diagram (Fig. 2). The approxi- 
mately hyperbolic shape of this plot suggests that to a 
first approximation, the source can be understood as 
the superposition of a constant and a varying source, 
each with nearly constant hardness ratio. An explicit 
example of this is seen in the data of Fig. 1 (bot- 
tom panel), where the hardness ratio varies in a near 
square-wave pattern as the count rate alternates be- 
tween a low and high value. The spectra in these two 
brightness levels are strikingly different (inset frame 
in Fig. 2), mainly because of a large excess of emis- 
sion below 20 keV at high count rate. This would 
appear to be the varying component invoked above. 
We confirmed this idea by the following procedure. 

We accumulated X-ray spectra for the quiescent 
and outburst intervals for the observation of 1996 Oc- 
tober 7th (Fig. 2, inset). Fits with simple single- 
component models did not give acceptable results. 
We obtained good fits (reduced 'x^=\.2 and 0.8 re- 
spectively) with a "standard" model for black-hole 
candidates, consisting of a multi-temperature disk - 
blackbody plus a power law ( Mitsuda et al. 1984 ). 
The addition of an iron line plus absorption edge was 
necessary to fit the data around 6 keV. All fits were 
performed in the range 2-50 keV, and a 1% system- 
atic uncertainty was added to the data to account for 
t he uncertainties in the response matrix calibration 
( Cui et al. 1997 ). The best fit parameters, using 
the known distance and inclination (we assume that 
the latter is the same as the inclination of the radio 
jets to the line of sight), are shown in Table 1. For 



the outburst state we find a column density Nh = 
8.0x10^2 cni-2, somewhat higher than the ROSAT 
value (6x10^^ cm~^, see Greiner, Morgan, & Remil- 
lard 1996). For the quiescent state we have fixed Nr 
to the outburst value. The marked spectral changes 
between quiescence and outburst are due to a steep- 
ening of the power law component and a large change 
in temperature of the disk-blackbody. From addi- 
tional fits to other parts of the light curve and to other 
observations we find that while the power law index 
can assume values anywhere in the range 2.3-3.9, the 
inner temperature of the disk-blackbody component 
switches between two values: ~0.5 keV, which cor- 
responds to an inner radius of several hundred kilo- 
meters, and ~2.2 keV, with an inner radius of ~20 
km. This hot disk component switches on at count 
rates larger than 20000 c/s (see Fig. 2), even during 
the brief peaks of the oscillating events in the flaring 
state. Below 20000 c/s, the hardness ratio changes in 
Fig. 2 are due primarily to changes in the power-law 
index. 

By looking at the light curve in Figure 1 at high 
time resolution, one sees that the fastest rise and fall 
times for the flare state remain constant. The shortest 
doubling rise times are ~2 s, while the corresponding 
shortest decay times are significantly less, i.e. 0.5 - 
1.0 s. Two further characteristic time scales apparent 
in the X-ray light curve (cf. Fig. 1) are the duration 
of the quiescent state and of a complete quicscent- 
outburst-flare cycle. For the light curve plotted in 
Fig. 1, these are ~300 s and ~ 1000 s, respectively. 

3. INTERPRETATION 

If we take the simplest possible interpretation of 
the multi-temperature blackbodies, i.e. as the emis- 
sion from optically thick regions of an accretion disk 
around a compact object of mass M, we have i?in, v — 
20 km for the brightest state of the varying source, 
and i?i„_ c ~300 km for the constant source. From 
R\n and T-^ we can infer the mass accretion rate ac- 
cording to 



M = 8TrRf^ aT ^j3GM. 



-01 



i?in,v cannot be smaller than the innermost stable or- 
bit around a black hole. This gives M < 2.4Mo, M > 
1 X IQ-'^ Mq yr^i in outburst and M > 6.8 x 10""^ M© 
yr~^ in quiescence for a Schwarzschild hole {Rin,v > 
GGM/c^). For an extreme Kerr hole we have Rin,v > 
GM/c^, so M < 14Mo, and M > 1.5 x 10"^ M© 



escence. The fact that the inferred accretion rates 
anticorrelate with the total luminosity implies that 
the accretion yield is lower when the accretion rate 
is high. This can be accomplished by trapping of the 
radiation and quasi-spherical infall ( "advection" ) of 
the accretion energy into a black hole (cf pK^atz 1977 , 
Begelman 197S, Narayan 1996|). We shall see that 



the disk interpretation indeed suggests that just this 
occurs. 

We arrive at the following picture. The "constant" 
source is the emission from the outer parts of the ac- 
cretion disk. This component varies only on the long 
viscous time scales characteristic of the outer disk. 
The varying source is the emission from the inner disk 
region and from a hot comptonizing corona above the 
disk, whose spectrum varies in response to changes in 
the disk continuum. As we shall explain below, the 
inner disk region empties and refills on time scales 
of seconds, thus causing the rapid flux and spectral 
changes. Radiation pressure becomes important in 
the inner region: it exceeds gas pressure for disk radii 



R< 10' 



8„0.1 



MiV^io 



0.14 



cm 



(2) 



(cf Frank, King, fc Raine 1992]). Here a is the 



Shakura-Sunyaev viscosity parameter, Mig is the ac- 
cretion rate in units of 10^^ g s^^, and Afio is the 
black hole mass in units of lOM©. In such regions the 
disk scale height H is almost constant with radius, so 
that 

^ ~ O.I5M18-R7' (3) 



(cf Frank, King, fc Raine 1992 ) and the thin disk 



approximation necessarily breaks down near the inner 
edge of the disk {R^ ~ 0.1). Here R^ is the radial 
coordinate in units of 10^ cm. 

Thermal balance at a given radius of a disk can 
be expressed in terms of a relation between the local 
accretion rate M and the surface density S. It is well 
known that for regions with Prad ~ ^gas, the M(S) 
curve has a branch with a negative slo pe, which is 
thus thermally and viscously unstable ( [Lightman & 
Eardley 1974 ). As pointed out by Abramowicz et 



yr 



"1 in outburst and Af > 1 x 10"'^ M, 



yr in qui- 



al. (1988), Lasota & Pelat (1991) and Chen et al. 
(1995), inward advection of the hot gas in the high 
M state is likely to stabilize the system. This leads 
to an S-shaped thermal balance curve in the Af — E 
plane, at typical disk effective temperatures ^ few 
X 10® K on the lower branch. This is very similar to 
the familiar thermal balance curve found near the hy- 
drogen ionization temperature ~ 10** K in dwarf nova 



disks. In the same way as there, the disk is locally 
thermally stable on the branches with positive slope, 
and locally unstable on the middle branch with nega- 
tive slope. We obviously again have the possibility of 
limit-cycle behavior, this time only in central regions 
of the disk with temperatures exceeding a few million 
degrees. Cannizzo (1996) has modeled the variability 
of the bursting pulsar GRO J1744-28 (Kouvcliotou 
et ; .l. 1996| ) in these terms. For that source the large 
assumed value of the inner disk radius (because of 
the neutron star magnetic field) leads to much longer 
time scales than for GRS 1915+105. We note that 
the Lightman-Eardley instability is only expected to 
operate inside the radius given by Eq. (^, which in- 
deed is close to the inner radius we estimate for the 
"constant" source. 

Since the upper branch of the S-curve involves 
physics which is currently unclear, we confine our- 
selves to estimating the time scales of the X-ray vari- 
ation. After an outburst, the outer disk attempts 
to refill the central disk via steady inflow on its vis- 
cous time scale, initially causing the slow recovery 
during the quiescent phase. A brightness rise occurs 
as a heating wave propagates outwards into the cen- 
tral disk from its inner edge before it has been refilled 
to a steady-state density profile. The hot radiating 
matter is removed by infall into the central black hole 
on the effective viscous time in the outbursting state. 
The latter time is not well known, as the thin disk 
approximation fails in this state. Heating fronts and 
infall can alternate many times in the flaring state, 
before the central disk eventually returns to the qui- 
escent state. The whole cycle repeats roughly on the 
viscous time of the outer disk. 

We can estimate the expected time scales as fol- 
lows. The rise time is the time for the heating wave 
to pass through the central disk. Assuming that the 
heating front travels with subsonic velocity, the sound 
speed cs being calculated in the cool (< 6 x 10^ K) 
disk state, we get 



time scale at the outer edge of the central disk, i.e. 



f ■ > 



Rn 



Cs 



l^out, 7 S. 



(4) 



where i?out, 7 = -Rout/10^ cm. The decay time cannot 
be shorter than the free-fall time scale, giving 



id( 



ccay 



.3/2 ..-1/2 
•-out, 7-"^ 10 



(5) 



The whole cycle should repeat roughly on the viscous 



_R^ 

^rec — 



R R 

H acs 



Qa-^M^^R^^ s. 



(6) 



where we have written the kinematic viscosity as 
V = aCgH and used (Eqs. and 0). We note that 
(Eqs. U and ^) are in rough agreement with the ob- 
served rise and decay times '^ 2 s, ~ 0.5 s, while 
typical values a ~ 10~^ assumed for quiescent disks 
imply recurrence times ^ QOOR'^ s, quite close to the 
observed values. Note that the X-ray rise and fall 
episodes differ from the case considered by Cannizzo 
(1996). There the infall on to the central neutron star 
is the process which raises the X-ray production rate, 
whereas here infall causes hot matter to be lost down 
the black hole, leading eventually to a decay of the 
X-rays. 

We note that the several features of the X-ray light 
curves strongly resemble the results of dwarf nova disk 
instability calculations. Even though the basic cause 
of the instability differs here, the general picture of 
heating and cooling fronts travelling back and forth 
leads us to expect generic similarities. In particular, 
alternating wide and narrow outbursts, and series of 
narrow outbursts followed by a wide one are com- 



mon features (see Van Paradijs 1983, for the outburst 
width distribution for dwarf novae) . These typically 
result when the heating wave stalls before crossing 
the entire unstable region. Similarly, a generic pre- 
diction of this type of model comes from the fact that 
the rise and decay episodes are fundamentally asym- 
metric. When observed with high enough time res- 
olution, the system should perform a hysteresis loop 
in an X-ray color-color diagram. Unfortunately, the 
large variations in the power law component do not 
currently allow one to draw a conclusion from sim- 
ple color-color plots. The accumulation of high time 
resolution energy spectra, not available in the current 
data, will allow one to test for the presence of such a 
hysteresis loop. A full quantitative discussion will be 
presented in a future paper. 

The extremely high values of the mass accretion 
rate we infer from the spectral fits in the quiescent 
state do not lead to a correspondingly high X-ray lu- 
minosity generation in the inner disk (with a canon- 
ical efficiency of 0.2 c^ per gram one would have ex- 
pected an X-ray luminosity in excess of 10^" erg/s). 
This indicates that a major fraction of the gravita- 
tional energy released in the inner disk is not emitted 
as radiation, suggestive of advection into a black hole. 



Assuming that the accretion rates derived from our 
spectral model are essentially correct and that ejec- 
tion is relatively unimportant, a substantial fraction 
of the inflow disappears into the central object with- 
out broadcasting a signal to the outside world. This 
supports the idea that it does not have a solid sur- 
face, and thus has the most characteristic property 
of a black hole. We note that the interpretation pro- 
posed in this paper depends on the spectral model 
we have used. We emphasize that at this point it is 
difficult to exclude absolutely the possibility that the 
compact object in GRS 1915+105 is a neutron star, as 
the study of X-ray spectra of neutron stars accreting 
at a hyper-Eddington rate is fairly unexplored. How- 
ever, a black hole surrounded by an optically thick ac- 
cretion disk of variable inner radius provides a natural 
explanation for the otherwise apparently intractable 
phenomenology presented by GRS 1915-1-105. 

It is improbable that the high inferred accretion 
rate represents the high state of some much longer- 
term instability in the outer disk. It seems more likely 
that the mass transfer rate is always very high, and 
that the outbursts of GRS 1915-1-105 are not caused 
by an increase of the mass transfer rate, but quite 
possibly by a decrease from a value high enough to 
smother the X-ray emission. In the smothered state, 
the object would then appear only as a bright infrared 
source. Average mass transfer rates of the required 
order (lO~^-lO~^M0/yr) can easily be supplied for 
example by a donor star crossing the Hcrtzsprung 
gap, and thus undergoing thermal-timescale radius 



expansion (Kolb et al. 1997) 
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Fig. 1.— Upper panel: 2.0-13.3 keV PCA light curve. 
Time zero corresponds to 1996 October 7th 6:36 UT. 
Lower panel: corresponding hardness ratio (13.3-60.0 
kev / 2.0-13.3 keV). Examples of the three states can 
be identified as follows. Quiescent state: T~400-800; 
outburst state: T~800-1200; flaring state: T~1200- 
1700. 

Fig. 2. — Hardness-intensity diagram for the data in 
Fig. 1. Rate and hardness ratio are defined as in 
Fig. 1. Inset: PCA energy spectra for the average 
quiescent (dots) and outburst (circles) states. 



Table 1 
Best-fit parameters for the quiescence and outburst spectra of October 7th. Nh is the 

COLUMN absorption IN 10^^ CM~^, KTin IS THE INNER TEMPERATURE OF THE DISK-BLACKBODY (iN KeV), 
Rin IS THE INNER RADIUS IN KM, T IS THE POWER-LAW PHOTON INDEX. THE (uNABSORBED) FLUXES ARE IN 

"^ IN THE 2-50 keV band. 



State 


Nh 


kT™ 


Rin 


r 


DBB Flux 


PL flux 


Quiescence 
Outburst 


8.0 
8.0±0.5 


0.57±0.03 
2.27±0.03 


319±9 
20.3±0.3 


2.22±0.05 
3.57±0.03 


{0.96±0.01)xl0-8 
{3.99±0.07)xl0-» 


(0.80±0.03)xl0-8 
{4.65±0.22)xl0-8 
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